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The present status of the diquark model for exclusive reactions at moderately 
large momentum transfer is reviewed. That model is a variant of the Brodsky- 
Lepage approach in which diquarks are considered as quasi-elementary constituents 
of baryons. Recent applications of the diquark model, relevant to high energy 
physics with electromagnetic probes, are discussed: electromagnetic form factors 
of baryons in both the space-like and the time-like region, photoproduction of 
mesons, two-photon annihilations into proton-antiproton pairs as well as real and 
virtual Compton scattering on which the main emphasis is laid. 



Exclusive processes at large momentum transfer are described in terms of 
hard scatterings among quarks and gluons p. In this so-called hard scatter- 
ing approach (HSA) a hadronic amplitude is represented by a convolution of 
process independent distribution amplitudes (DA) with hard scattering ampli- 
tudes to be calculated within perturbative QCD. The DAs specify the distri- 
bution of the longitudinal momentum fractions the constituents carry. They 
represent Fock state wave functions integrated over transverse momenta. The 
convolution manifestly factorizes long (DAs) and short distance physics (hard 
scattering). The HSA has two characteristic properties, the power laws and 
the helicity sum rule. The first property says that, at large momentum transfer 
and large Mandelstam s, the fixed angle cross section of a reaction AB — > CD 
behaves as 

da/dt = f(9) s 2 - n (1) 
where n is the minimum number of external particles in the hard scattering 
amplitude. The laws (^) are modified by powers of logs. They also apply to 
form factors: a baryon form factor behaves as 1/Q 4 , a meson form factor as 
1/Q 2 . The counting rules are found to be in surprisingly good agreement with 
experimental data. Even at momentum transfers as low as 2 GeV the data 
seem to respect the counting rules. 

The second characteristic property of the HSA is the conservation of hadronic 
helicity. For a two-body process the helicity sum rule reads 

Xa + ^b = Ac + Ap. (2) 
It appears as a consequence of utilizing the collinear approximation and of 
dealing with (almost) massless quarks which conserve their helicities when in- 
teracting with gluons. The collinear approximation implies that the relative 
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Figure 1: The Pauli form factor of the proton scaled by Q°. Rata are taken from El. The 
solid line represents the result obtained with the diquark modelf]. 

orbital angular momentum between the constituents has a zero component in 
the direction of the parent hadron. Hence the helicities of the constituents 
sum up to the helicity of their parent hadron. The helicity sum rule is violated 
by 20 — 30% by many experimental data. A particular striking example is the 
Pauli form factor of the proton which is measured to be large □. Its Q 2 depen- 
dence (see Fig. 1) is compatible with a higher twist contribution (~ 1/Q 6 ). 
In explicit applications of the HSA (carried through only in leading twist and to 
lowest order QCD with very few exceptions) one encounters the difficulty that 
the data are available only at moderately large momentum transfer, a region 
in which non-perturbative dynamics may still play a crucial role. A general 
feature of such applications is the extreme sensitivity to the DAs chosen for 
the involved hadrons. Only strongly end-point concentrated DAs provide re- 
sults which are at least lor the magnetic form factor of the nucleon in fair 
agreement with the datad. This apparent success of the HSA is only achieved 
at the expense of strong contribution from soft regions where one of the con- 
stituents carries only a tiny fraction of its parent hadron's momentum. This 
is a very problematical situation for a perturbative calculation. It should be 
stressed that none of the DAs used in actual applications leads to a successful 
description of all large momentum transfer processes investigated so far. 
It seems clear from the above remarks that the HSA at leading twist although 
likely to be the correct asymptotic picture for exclusive reactions, needs rjiocL 
ifications at moderately large momentum transfer. In a series of papers □ l£j 
such a modification has been proposed by us in which baryons are viewed as 
composed of quarks and diquarks. The latter are treated as quasi-elementary 
constituents which partly survive medium hard collisions. Diquarks are an 
effective description of correlations in the wave functions and constitute a par- 
ticular model for non-perturbative effects. The diquark model may be viewed 
as a variant of the HSA appropriate for moderately large momentum transfer 
and it is designed in such a way that it evolves into the standard pure quark 



HSA asymptotically. In so far the standard HSA and the diquark model do 
not oppose each other, they are not alternatives but rather complements. The 
existence of diquarks is a hypothesis. However, from experimental and theo- 
retical approaches there have been many indications suggesting the presence 
of diquarks. For instance, they were introduced in baryon spectroscopy, in 
nuclear physics, in astrophysics, in jet fragmentation and in weak interactions 
to explain the famous AI = 1/2 rule. Diquarks also provide a natural expla- 
nation of the equal slopes of meson and baryon Regge trajectories. For more 
details and for references, seed. It is important to note that QCD provides 
some attraction between two quarks in a colour {3} state at short distances as 
is to be seen from the static reduction of the one-gluon exchange term. 
Even more important for our aim, diquarks have also been found to play a role 
in inclusive hard scattering reactions. The most obvious place to signal their 
presence is deep inelastic lepton-nuclean scattering. Indeed the higher twist 
contributions, convincingly observed LU, can be modelled as lepton-diquark 
elastic scattering. Baryon production in inclusive pp collisions-also reveals the 
need for diquarks scattered elastically in the hard interactiont3. For instance, 
kinematical dependences or the excess of the proton yield over the antiproton 
yield find simple explanations in the diquark model. No other explanation of 
these phenomena is known as yet. 

The diquark model: As in the standard HSA a helicity amplitude for the re- 
action AB — > CD is expressed as a convolution of DAs and hard scattering 
amplitudes (s, — t, —u ^> mf) 

M(s,t) = J &x c &x D dx A <lxB<S>c{xcW D {x D )T H {x u s,t)§ A {x A )§B{xB) (3) 

where helicity labels are omitted for convenience. Implicitly it is assumed in 
(||) that the valence Fock states consist of only two constituents, a quark and a 
diquarks (antiquark) in the case of baryons (mesons). In so far the specification 
of the quark momentum fraction Xi suffices; the diquark (antiquark) carries 
the momentum fraction 1 — Xi. If an external particle is point-like, e. g. a 
photon, the accompanying DA is to be replaced by 8(1 — Xi). Because of 
QCD evolution the DAs depend logarithmically on the momentum transfer. 
This fact is of minor importance in the limited range of momentum transfer 
in which data are available and is therefore ignored. As in the standard HSA 
contributions from higher Fock states are neglected. This is justified by the 
fact that that such contributions are suppressed by powers of a s /t as compared 
to those from the valence Fock state. 

In the diquark model spin (S) and spin 1 (V) colour antitriplet diquarks are 
considered. Within flavour SU(3) the S diquarks form an antitriplet, the V 
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diquarks an sixtet. Assuming zero relative orbital angular momentum between 
quark and diquark and taking advantage of the collinear approximation, the 
valence Fock state of an ground state octet baryon B with helicity A and 
momentum p can be written in a covariant fashion (omitting colour indices) 

\B;p,X) = fs^(x)B s u(p,\) + f v ^(x)B v ( 1 a + P a /m B ) To u(p,\)/V3 (4) 

where u is the baryon's spinor. The two terms in (Q) represent configurations 
consisting of a quark and either a scalar or a vector diquark, respectively. The 
couplings of the diquarks with the quarks in a baryon lead to flavour functions 
which e. g. for the proton read 

B s = uS [u , d ] B v = [uV{ u4 ] - V%d V{ UiU y]/V3 ■ (5) 

The DAs &s(v) are conventionally normalized as J dxQ = 1. The constants 
fs(v) play the role the configuration space wave function at the origin. 
The DAs containing the complicated non-perturbative bound state physics, 
cannot reliably be calculated from QCD at present. It is still necessary to 
parameterize the DAs and to fit the eventual free parameters to experimental 
data. Hence, both the models, the standard HSA as well as the diquark model, 
only get a predictive power when a number of reactions involving the same 
hadrons are investigated. In the diquark model the following DAs have been 
proven to work satisfactorily well in many applicationsQ 

$f (x)=Nix(l - x) 3 exp [-b 2 (m 2 q /x + m|/(l - x))] (6) 
$£(x) =Nyx{l - xf{l + 5.8 x - 12.5 x 2 ) exp [~b 2 (m 2 q /x + m 2 v /(l - x))] . 

These DAs are a suitable adaption of a meson DA obtained by transforming 
the harmonic oscillator wave function to the light-cone. The constants N B are 
fixed through the normalization convention (e. g. for the proton Ng = 25.97 
and Ny = 22.92). The DAs exhibit a mild flavour dependence via the expo- 
nential which also guarantees a strong suppression of the end-point regions. 
The masses in (|^) are constituent masses since they enter through a rest frame 
wave function. For u and d quarks we take 350 MeV and for the diquarks 
580 MeV. Strange quarks and diquarks are assumed to be 150 MeV heavier 
that the non-strange ones. It is to be stressed that the quark and diquark 
masses only appear in the DAs (]|); in the hard scattering kinematics they are 
neglected. The final results (form factors, amplitudes) depend on the actual 
mass values mildly. The transverse size parameter b is fixed from the assump- 
tion of a Gaussian transverse momentum dependence of the full wave function 
and the requirement of a value of 600 MeV for the mean transverse momen- 
tum (actually b = 0.498 GeV -1 ). As the constituent masses the transverse 
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size parameter is not considered as a free parameter since the final results only 
depend on it weakly. 

The hard scattering amplitudes Tjj determined by short-distance physics, are 
calculated from a set of Feyman graphs relevant to a given process. Diquark- 
gluon and diquark-photon vertices appear in these graphs which, following 
standard prescriptions, are defined as 

SgS : ig s t a (p 1 +p 2 )^ 

VgV : -i g s t a {g a p{pi + P2)» ~ 9f)^ [(1 + «) P2 ~ K Pi] a 

~9na [(1 + k) pi - np 2 ] p | (7) 

where g s = y 'Ana s is the QCD coupling constant, n is the anomalous magnetic 
moment of the vector diquark and t a = X a /2 the Gell-Mann colour matrix. For 
the coupling of photons to diquarks one has to replace g s t a by — y/Awaen where 
a is the fine structure constant and eu is the electrical charge of the diquark 
in units of the elementary charge. The couplings DgD are supplemented by 
appropriate contact terms required by gauge invariance. 

The composite nature of the diquarks is taken into account by phenomenolog- 
ical vertex functions. Advice for the parameterization of the 3-point functions 
(diquark form factors) is obtained from the requirement that asymptotically 
the diquark model evolves into the standard HSA. Interpolating smoothly be- 
tween the required asymptotic behaviour and the conventional value of 1 at 
Q 2 = 0, the diquark form factors are actually parametrized as 

The asymptotic behaviour of the diquark form factors and ±h.e connection to 
the hard scattering model is discussed in more detail in Ref. Em. In accordance 
with the required asymptotic behaviour the n-point functions for n > 4 are 
parametrized as 

av WT&) F ^ m - (9) 

The constants as,v are strength parameters. Indeed, since the diquarks in 
intermediate states are rather far off-shell one has to consider the possibility 
of diquark excitation and break-up. Both these possibilities would likely lead 
to inelastic reactions. Therefore, we have not to consider these possibilities ex- 
plicitly in our approach but excitation and break-up lead to a certain amount 
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of absorption which is taken into account by the strength parameters. Admit- 
tedly, that recipe is a rather crude approximation for n > 4. Since in most 
cases the contributions from the n-point functions for n > 4 only provide small 
corrections to the final results that recipe is sufficiently accurate. 
Special features of the diquark model: The diquark hypothesis has striking con- 
sequences. It reduces the effective number of constituents inside baryons and, 
hence, alters the power laws (0). In elastic baryon-baryon scattering, for in- 
stance, the usual power s~ 10 becomes s~ 6 F(s) where F represents the effect 
of diquark form factors. Asymptotically F provides the missing four powers of 
s. In the kinematical region in which the diquark model can be applied (—t, 
—u> 4GeV 2 ), the diquark form factors are already active, i. e. they supply a 
substantial s dependence and, hence, the effective power of s lies somewhere 
between 6 and 10. The hadronic hclicity is not conserved in the diquark model 
at finite momentum transfer since vector diquarks can flip their helicities when 
interacting with gluons. Thus, in contrast to the standard HSA spin-flip de- 
pendent quantities like the Pauli form factor of the nucleon can be calculated. 
Electromagnetic nucleon form factors: This is the simplest application of the 
diquark model and the most obvious place to fix the various parameters of the 
model. The Dirac and Pauli form factors of the nucleon are evaluated from 
the convolution formula (||) with the DAs (Q) and the parameters are deter- 
mined from a best fit to the data in the space-like region. The following set of 
parameters 

fs = 73.85 MeV, Q| = 3.22GeV 2 , a s = 0.15, ( ] 

fv = 127.7 MeV, Q\ = 1.50 GeV 2 , a v = 0.05, k = 1.39 ; ' 

provides a good fit of the dataB. a s is evaluated with Kqcd — 200 MeV and 
restricted to be smaller than 0.5. The parameters Qs and Qy, controlling the 
size of the diquarks, are in agreement with the higher-twist effects—observed 
in the structure functions of deep inelastic lepton-hadron scattering til if these 
effects are modelled as lepton-diquark elastic scattering. The Dirac form factor 
of the proton is perfectly reproduced. The results for the Pauli form factor are 
shown in Fig. 1. The predictions for the two neutron form factors are also in 
agreement with the data. However, more accurate neutron data are needed in 
the Q 2 region of interest in order to examine the model crucially. The nucleon's 

r:ial form factorQ and its electromagnetic form factors in the time-like regions 
have also been evaluated. Both the results compare well wilL data. Even 
electroexcitation of nucleon resonances has been investigated!!^!!]. 
Real Compton scattering (RCS): jp — > jp is the next reaction to which the 
diquark model is applied. Since again the only hadrons involved are protons 
RCS can be predicted in the diquark model; there is no free parameter to 
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Figure 2: Typical Feynman graphs contributing to yi*' p — » 7 p. 
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Figure 3: The cross section for RCS off protons scaleci— by s 6 vs. cos(9 for three different 
photon energies. The experimental data are taken from I 1 ^. 



be adjusted. Typical Feynman graphs contributing to that process are shown 
in Fig. 2. The results of the diquaxk model for RCS are shown in Fig. 3 
for three different photon energies tffi. Note that in the very forward and 
backward regions the transverse momentum of the outgoing photon is small 
and, hence, the diquark model which is based on perturbatpe QCD, is not 
applicable. Despite the rather small energies at which dataEj are available, 
the diquark model is seen to work rather well. The predicted cross section 
does not strictly scale withes -6 . The results obtained within the standard 
HSA are of similar quality E£l. The diquark model also predicts interesting 
photon asymmetries and spin correlation parameters (see the discussipji inQ). 
Even a polarization of the proton, of the order of 10%, is predicted Q. This 
comes about as a consequence of helicity flips generated by vector diquarks 
and of perturbative phases produced by propagator poles appearing within 
the domains of momentum fraction integrations. The poles are handled in 
the usual way by the is presription. The appearance of imaginary parts to 
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Figure 4: The integrated my — » pp cross section (| cos 6 |>-C.6). The solid line represents the 
diquark model prediction cl Data are taken from CLEOli3. 

leading order of a s is a non-trivial prediction of perturbative QCD 0; it is 
characteristic of the HSA and is not a consequence of the diquark hypothesis. 
Two-photon annihilation into pp pairs: This process is related to RCS by 
crossing, i. e. the same set of Feynman graphs contributes (see Fig. 2). The 
only difference is that now the diquark form factors are needed in the time-like 
region. The expressions (@,^) represent an effective parameterization of them 
valid at large space-like Q . Since the exact dynamics of the diquark system is 
not known it is not possible to continue these parameterizations to the time- 
like region in a unique way. A continuation can be defined as follows u: Q 2 is 
replaced by — s in (^,^) guaranteeing the correct asymptotic behaviour and, in 
order to avoid the appearance of unphysical poles at low Q 2 , the diquark form 
factors are kept constant once their absolute values have reached Cq = 1.3 u. 
The same definition of the time- like diquark form factors is used in the analysis 
the proton form factor in the time-like region. The diquark model predictions 
for the integrated 77 — > pp cross section is compared to the CLEO datalia in 
Fig. 4. At large energies the agreement between predictions and experiment is 
good. The predictions for the angular distributions are in agreement with the 
CLEO data too. The standard HSA on the other hand predicts a cross section 
which lies about an order of magnitude below the data El Recently CLEO 
has also measured two-photon annihilations into AA pairsEj. Surprisingly the 
integrated cross section is, within errors, as large as that for annihilations 
into pp pairs. Using the SU(6)-like spin-flavour dependence (||^), the diquark 
model predicts a AA cross section which is about a factor of 2 smaller than the 
CLEO data. The reason for this discrepancy is not yet understood. 
Virtual Compton scattering (VCS): This process is accessible through ep — > 
epj. An interesting element in that reaction is that, besides VCS, there is 
also a contribution from the Bethe-Heitler (BH) process where the final state 
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Figure 5: The cross section for VCS vs. cosd for several values of Q 2 /s at s = 5 GeV 2 . Upper 
left: the transverse cross section scaled by s 6 . Upper right: the ratio of the longitudinal over 
the transverse cross sections. Lower left (right): the ratio of the longitudinal (transverse) - 
transverse interference term over the transverse cross section. 

photon is emitted from the electron. Electroproduction of photons offers many 
possibilities to test details of the dynamics: One may measure the s, t and Q 2 
dependence as well as that on the angle 4> between the hadronic and leptonic 
scattering planes. This allows to isolate cross sections for longitudinal and 
transverse virtual photons. One may also use polarized beams and targets and 
last but not least one may measure the interference between the BH and the 
VC contributions. The interference is sensitive to phase differences. 
At s, —t and — u ^S> m 2 (or small |cos#| where 8 is the scattering angle of 
the outgoing photon in the photon-proton center of mass frame) the diquark 
model can also be applied to VCS □. Again there is no free parameter in 
that calculation. The relevant Feynman graphs are the same as for RCS (see 
Fig. 2). The model can safely be applied for s > 10 GeV 2 and | cos#| < 0.6. 
For the future CEBAF beam energy of 6 GeV the model is at its limits of 
applicability. However, since the diquark model predictions for real Compton 
scattering do rather well agree with the data even at s > 5 GeV 2 (see Fig. 3) 
one may expect similarly good agreement for VCS. Predictions for the VCS 
cross section are shown in Fig. 5. The transverse cross section (which, at 
Q 2 = 0, is the cross section for RCS is the dominant piece. The other cross 
sections only become sizeable for large values of | cos#|. Examination of the 
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Figure 6: The electron asymmetry in ep — > ep7 as predicted by the diquark modelEl. 

Bethe-Heitler contribution to the process ep — > epr/ reveals that it is small as 
compared to the VCS contribution at high energies, small values of | cos6>| and 
for an out-of-plane experiment, i. e. <f> > 50°. 

The last observable I want to discuss is the electron asymmetry in ep — > epj: 



A L = 



g(+)-gR 
a(+)+a(+) 



(11) 



where ± indicates the helicity of the incoming electron. Al measures the 
imaginary part of the longitudinal - transverse interference. The longitudinal 
amplitudes for VCS turn out to be small in the diquark model (hence A\ c is 
small). However, according to the model, Al is large in the region of strong 
BH contamination (see Fig. 6). In that region, A^ measures the relative phase 
(being o£_Derturbative origin from on-shell going internal gluons, quarks and 
diquarksliH) between the BH amplitudes and the VCS ones. The magnitude 
of the effect shown in Fig. 6 is sensitive to details of the model and, therefore, 
should not be taken literally. Despite of this our results may be taken as an 
example of what may happen. The measurement of Al, e. g. at CEBAF, will 
elucidate strikingly the underlying dynamics of VCS. 

Photo- and electroproduction of mesons: This is already a quite complicated 
reaction to which all together 158 Feynman graphs contribute._|Jp to now only 
the two processes jp — > K + A, K* + A have been investigated EH. The analyses 
of other final states as well as electroproduction are in progress. The calcu- 
lation of KA production is somewhat simpler than that for other final states 
because only scalar diquarks contribute. The analysis of many different final 
states will provide deep insight in the dynamics. 

As compared to the processes discussed above a new element_appears now, 
namely the mesonic DA. Comparison of predictions with data til (see Fig. 7) 
revealed that the asymptotic form of the Kaon DA (~ x(l — x)) works very well 
(using the standard value of the Kaon decay constant). On the other hand the 
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Figure 7: Differential cross section for 7fl—r+ A" "'"A scaled by s 7 vs. cos6 cm . Solid (dash- 
double-dotted) line: Diquark model result|]J at p^ ab = 6(4) GeV using the asymptotic Kaon 
DA. Short-dashed line: Results _.= 6 GeV using the CZ Kaon DA. Long-dashed line: 
Predictions from the standard HSAEj. The data are taken from Anderson et al. Ej- 

double-humped DA proposed by Chernyak and Zhitnitski! fails, the predicted 
cross section is ±00 large as compared with the data. The predictions from the 
standard HSA E3 are smaller than those from the diquark model. It should 
be mentioned that for the A the SU(6)-like spin-flavour dependence ||) is 
used. How to reconcile this with the apparent failure of the SU(6)-like A wave 
function in 77 — ► AA remains to be seen. 

Summary and outlook: The diquark model which represents a variant of the 
HSA, combines perturbative QCD with non-perturbative elements. The di- 
quarks represent quark-quark correlations in baryon wave functions which 
are modelled as quasi-elementary constituents. This model has been applied 
to many photon induced exclusive processes at moderarely large momentum 
transfer (typically > 4 GeV 2 ). From the analysis of the nucleon form factors the 
parameters specifying the diquark and the DAs, are fixed. Compton scatter- 
ing and two- photon annihilations of pp can then be predicted. The comparison 
with existing data reveals that the diquark model works quite well and in fact 
much better then the pure quark HSA. Using the asymptotic DA for the Kaon 
and SU(6) ideas to fix the A DA one can also predict photoproduction of KA. 
Again there is agreement between predictions and experiment. 
Predictions for the VCS cross section and for the ep — > epj cross section have 
also been made for kinematical situations accessible at the upgraded CEBAF 
and perhaps at future high energy accelerators like ELFEQHERA. According 
to the diquark model the BH contamination of the photon electroproduction 
becomes sizeable for small azimuthal angles. The BH contribution also offers 
the interesting possibility of measuring the relative phases between the VC 
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and the BH amplitudes. The phases of the VC amplitudes are a non-trivial 
phenomenon generated by the fact that some of the internal quarks, diquarks 
and gluons may go on mass shell. The electron asymmetry Al is particularly 
sensitive to relative phases. In contrast to the standard HSA the diquark model 
allows to calculate helicity flip amplitudes, the helicity sum rule (|) does not 
hold at finite Q 2 . One example of an observable controlled by helicity flip con- 
tributions is the Pauli form factor of the proton. Also in this case the diquark 
model accounts for the data. 



References 

1. G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

2. P. Bosted et al., Phys. Rev. Lett. 68, 3841 (1992). 

3. V. L. Chernyak and A. R. Zhitnitsky, Phys. Rep. 112, 173 (1984). 

4. M. Anselmino, P. Kroll and B. Pire, Z. Phys. C 36, 89 (1987). 

5. P. Kroll, Proceedings of the Adriatico Research Conference on Spin and 
Polarization Dynamics in Nuclear and Particle Physics, Trieste, 1988. 

6. P. Kroll, W. Schweiger and M. Schurmann, Int. Jour, of Mod. Physics 
A 6, 4107 (1991) 

7. R. Jakob, P. Kroll, M. Schurmann and W. Schweiger, Z. Phys. A 347, 
109 (1993). 

8. P. Kroll, Th. Pilsner, M. Schurmann and W. Schweiger, Phys. Lett. B 
316, 546 (1993). 

9. P. Kroll, M. Schurmann and P. Guichon, Nucl. Phys. A 598, 435 (1996). 

10. P. Kroll, M. Schurmann, K. Passek and W. Schweiger, preprint UNI- 
GRAZ - UTP 15-04-96. 

11. M. Virchaux and A. Milsztajn, Phys. Lett. B 274, 221 (1992). 

12. A. Breakstone et al, Z. Phys. C 28, 335 (1985). 

13. P. Kroll, M. Schurmann and W. Schweiger, Z. Phys. A 342, 429 (1992). 

14. J. Bolz, P. Kroll and J. G. Korner, Z. Phys. A 350, 145 (1994). 

15. M. A. Shupe et al., Phys. Rev. D 19, 1921 (1979). 

16. A. S. Kronfeld and B. Nizic, Phys. Rev. D 44, 3445 (1991). 

17. G. R. Farrar et al., Phys. Rev. Lett. 62, 2229 (1989). 

18. M. Artuso et al., CLEO collaboration, Phys. Rev. D 50, 5484 (1994). 

19. G. R. Farrar, E. Maina and F. Neri, Nucl. Phys. B 259, 702 (1985); B 
263, 746 (1986)(E). 

20. D. W. Bliss, representing the CLEO collaboration, APS meeting 1996. 

21. R. L. Anderson et al., Phys. Rev. D 14, 679 (1976). 

22. G. R. Farrar, K. Huleihel and H. Zhang, Nucl. Phys. B 349, 655 (1991). 



12 



